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Winding number transitions in the two dimensional softly broken 0(3) nonlinear sigma 
model are studied at finite energy and temperature. New periodic instanton solutions which 
dominate the semiclassical transition amplitudes are found analytically at low energies, and 
numerically for all energies up to the sphaleron scale. The Euclidean period /? of these 
finite energy instantons increases with energy, contrary to the behavior found in the abelian 
Higgs model or simple one dimensional systems. This results in a sharp crossover from 
instanton dominated tunneling to sphaleron dominated thermal activation at a certain critical 
temperature. Since this behavior is traceable to the soft breaking of conformal invariance 
by the mass term in the sigma model, semiclassical winding number transition amplitudes 
in the electroweak theory in 3 + 1 dimensions should exhibit a similar sharp crossover. We 
argue that this is indeed the case in the standard model for Mh < 4:Mw- 
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I. INTRODUCTION 

Gauge theories of the strong and electroweak interac- 
tions are characterized by a multiple vacuum structure. 
Tunneling transitions between different vacua are respon- 
sible for physically interesting effects, such as baryon 
number violation in the electroweak theory At zero 
temperature and energy these winding number transi- 
tions are dominated by the familiar zero energy instanton 
solutions of the Euclidean field equations, with vacuum 
boundary conditions ||^. 

At finite temperatures thermal activation over the po- 
tential barrier separating the multiple vacua can take 
place, in addition to quantum tunneling. The static clas- 
sical solution whose energy is equal to the top of this bar- 
rier between neighboring vacua is the sphaleron j^]. At 
sufficiently high temperatures, transitions between differ- 
ent winding number sectors are dominated not by quan- 
tum tunneling but by classical thermal activation, with 
a rate controlled by the energy of the sphaleron |^ . 

In simple one dimensional quantum mechanics, as the 
temperature is increased there is typically a smooth 
crossover from zero energy quantum tunneling via the 
instanton to high temperature thermal activation via 
the sphaleron Q. The corresponding classical solutions 
which interpolate between these two situations are known 
as periodic instantons (or "calorons") These finite 

energy solutions of the classical Euler-Lagrange equations 
possess turning points at finite Euclidean time P, and 
dominate the semiclassical transition rate at a temper- 
ature related to the Euclidean period by ksT = h/f3. 
(Hereafter we set Ti — — c — 1, so that T — f3^^.) 

Similar considerations are expected to apply to finite 
energy winding number transitions in quantum field the- 



ory, although there the situation is much less well ex- 
plored and very few of the required Euclidean periodic 
classical solutions are known. In particular, it is not clear 
a priori if the smooth crossover of the transition rate from 
quantum tunneling to thermal activation is a generic fea- 
ture, or if not, on what aspects of the field theory this 
behavior depends. Our purpose in this paper is to inves- 
tigate this question by finding and studying the finite en- 
ergy periodic instanton solutions in a specific model, the 
0(3) nonlinear sigma model in two dimensions, modified 
by a suitable conformal symmetry breaking term. This 
model is chosen for the many features it shares in com- 
mon with spontaneously broken gauge theories in four 
dimensions, and in particular, the standard model of elec- 
troweak interactions. In the 5/7(2) x U{1) electroweak 
theory the rate of baryon number violating winding tran- 
sitions at high energies remains an open question, despite 
considerable efforts in recent years A detailed study 
of the solutions of the classical nonlinear field equations 
of the theory appears to offer the best hope of addressing 
this issue ]9|,p^. 

The 0(3) sigma model possesses the advantage of be- 
ing simple enough that the zero energy instanton, the 
finite energy sphaleron, and the spectrum of linearized 
fluctuations about each of these solutions are all known 
analytically. These analytic results are fixed markers in 
the classical solution space that can be used as launching 
points for a detailed numerical study of the finite energy 
periodic instanton solutions, and which provide useful 
non-trivial checks on the numerical methods. The nu- 
merical techniques developed in the context of the sigma 
model may be applied then with considerably more con- 
fidence to four dimensional gauge theories such as the 
electroweak theory. The present work may be regarded 
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as the first step in this program. 

The chief result of our study of periodic instanton so- 
lutions in the softly broken 0(3) sigma model is that the 
Euclidean period (3 of the finite energy instanton solu- 
tions increases with increasing energy, i.e. 



d(3{E) 
dE 



> 0, 



(1.1) 



for all energies from zero up to the sphaleron. Since 
the Euchdean solution with period /? is associated with 
the winding number transition amplitude at temperature 
T = /?^^, one might expect that it should be T which in- 
creases with energy, and not (3. Indeed in simpler models 
such as the one dimensional quantum pendulum or the 
two dimensio nal abelian Higgs model the very opposite 
behavior to (1.1) is found. Our analytic and numerical 
study of perio dic i nstanton solutions in the 0(3) model 
will show that (1.1) holds all the way up to the sphaleron 
energy, where the periodic instantons merge with the 
sphaleron. Continuing further to energies E > Egpji re- 
sults in the periodic instanton solutions moving off into 
the complex domain with still increasing real period. 

As we shall see, the main physical consequence of this 
increase of period with energy is that there is a sharp 
crossover from instanton dominated quantum tunneling 
to sphaleron dominated thermal activation at a certain 
critical temperature (of the order of the mass in the 
model), rather than a smooth transition between the 
two. Since this feature is traceable to the conformal in- 
variance of the unbroken 0(3) model, the same sharp 
crossover from quantum tunneling to thermal activation 
should be expected in four dimensional gauge theories 
where the conformal and spontaneous gauge symmetry 
breaking arises from the Higgs sector vacuum expecta- 
tion value. Thus, the global picture which emerges from 
our study of the space of classical Euclidean solutions 
and their contribution to fixed energy and fixed tem- 
perature winding number transitions in the 0(3) model 
should be directly applicable to the electroweak theory. 
We present evidence that this is indeed the case, at least 
for not too large Higgs' self-coupling, i.e. for Higgs' mass 
Mh < ^Mw 

The structure of the paper is as follows. In the next 
section we discuss the general properties of periodic in- 
stanton solutions and their contribution to winding num- 
ber transition rates at finite energy and temperature, in 
the context of a quantum mechanical model with only one 
degree of freedom, viz. the simple pendulum. In Section 
3 we review the 0(3) nonlinear sigma model both before 
and after introducing the explicit soft symmetry break- 
ing mass term. The mass term is necessary in order for 
a finite energy sphaleron solution to exist. In Section 4 
we consider finite energy instantons in the softly broken 
sigma model, employing analytic techniques at low en- 
ergies where perturbation theory about the zero energy 
instantons of the unbroken model is applicable, and dis- 
cussing the behavior of the periodic solutions near the 



sphaleron energy E^ph- In Section 5 we describe our nu- 
merical techniques and present results for the finite en- 
ergy instanton solutions smoothly interpolating between 
the two limits. We conclude in Section 6 with a discussion 
of our results, and their implications for the analogous 
winding number transitions in the electroweak theory. 



II. PERIODIC INSTANTONS IN QUANTUM 
MECHANICS 

At finite temperature the partition function for a quan- 
tum system may be expressed as a path integral over 
configurations with fixed Euclidean periodicity 



Z(/3) 



9(0)=9(/3) 



[Vq] exp{-S[q;(3]) 



(2.1) 



where S is the classical Euclidean action evaluated on 
the configuration over one period p. If there is a small 
coupling in the problem which can be scaled out of the 
action hy S = s/g^ and we consider the weak coupling 
limit —^ 0, then the path integral is dominated by 
the extrema of the classical action, i.e. solutions of the 
classical Euclidean equations of motion obeying periodic 
boundary conditions. 

Strictly speaking only stable solutions, i.e. local min- 
ima of the action can contribute to the true equilibrium 
partition function at finite temperature. However, we are 
interested here in calculating the rate of real time wind- 
ing number transitions at finite temperature or energy. 
In order to contribute to such a transition rate a classical 
solution should not be a strict minimum but rather a sad- 
dle point possessing exactly one negative mode direction. 
Analytic continuation of the semiclassical approximation 
to Z{f3) in this one negative mode direction leads to an 
imaginary part which may be interpreted as a real time 
transition rate fll iP . A similar conclusion follows from an 



analytic extension of the path integral in (2.1) to complex 
valued configurations |1^ . 

In systems with only one degree of freedom it is 
straightforward to find the relevant periodic instanton 
solutions with a single negative mode direction by sim- 
ple quadrature. Since the energy is an integral of the 
motion we have simply to calculate 



PiE) 



dq 



^2V{q) - 2E 



(2.2) 



over the periodic trajectory with fixed energy E begin- 
ning and ending at the same turning point of the poten- 
tial V{q). The Euclidean action corresponding to this 
periodic trajectory is 



S{P)^jdq V2V{q)-2E{f3) 



(2.3) 



where E{(3) is obtained by inverting (|2.2|) 
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In order to make the discussion definite let us consider 
the case of a simple pendulum with the classical periodic 
potential, 



V{q) — U)'^ {1 — cos q ) 



(2.4) 



The periodic local minima of this potential at q = 2mr 
correspond to the periodic ground state vacuum structure 
we find in non-abelian gauge theories. The local maxima 
at g = (2n + 1)tt are the static but unstable "sphaleron" 
solutions of this simple model with energy, 



Esph — 2w 



(2.5) 



In this case the ratio of the zero point energy, tj/2 of 
the harmonic potential in the vicinity of the minima to 
the height of the potential barrier between minima, i.e. 
Lo/2Esph = (4w)~^ is the parameter which must be small 
(in units of ?i = / = 1) in order to justify a semiclassi- 
cal treatment of the quantum pendulum with moment of 
inertia / about its pivot point. 

At very low temperatures and energies the winding 
number transitions between neighboring minima at and 
27r are dominated by quantum tunneling with a rate of 
order exp(— 26*0) where Sq is the classical action for the 
zero energy kink solution. 



cos 



<?o(t) 



' tanh(a;T) 



(2.6) 



which interpolates between the two minima in infinite 
Euclidean time r. The anti-kink solution is obtained from 
(2.6) by time reversal r — ^ — r and the periodic trajectory 
beginning and ending at the same vacuum is a widely 
separated kink/anti-kink pair with total action 25*0 = 
16u; 3> 1, which implies that tunneling is exponentially 
suppressed at zero temperature or energy. 

At the same time it is quite clear that this exponen- 
tial suppression disappears at temperatures or energies 
comparable to Esph when the pendulum can jump over 
the barrier between the neighboring minima by classi- 
cal thermal activation. The semiclassical behavior of the 
transition rate at intermediate energies and temperatures 
is easily found analytically i n th e case of t he pendulum 
since both of the integrals (2.2) and (2.3) may be ex- 
pressed in terms of the complete elliptic integrals K and 
E as 



/3= -K(fc) 

UJ 



and 



5 = 8w [2E(fc) - (1 - fc2)K(fc)] , 



(2.7) 



(2.8) 



where k is the modulus of the elliptic functions, related 
to the energy by 



or 



k = 



The actual periodic instanton solution q{T) with these 
properties may be expressed in terms of a Jacobian ellip- 
tic function, sn, with modulus k by 



cos I ^^-r^ 1 = — fc sn(a;T, k) 



-tanh(tJT), E^Q 
0, E^ 2uP- 



(2.10) 



The corresponding behaviors of the Euclidean period and 
action in these limits are 



£:(/?) = 2u?k!' 



0, cx) 
Esphi (3 ^ P- 



and 



16(jj, /3 00 



where 



2^ 



(2.11) 



(2.12) 



(2.13) 



and k' = y/\ — k"^ is the complementary modulus of the 
elliptic functions. 
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FIG. 1. Exact action, Eqn . (2.8), (solid curve) and per- 
turbative action, Eqn. (2.26), (dashed curve) of the peri- 
odic instanton solutions of the simple pendulum model as a 
function of their Euclidean period {3. The curves start at 
/3 = /3_ = 27r/c^. 

In the first limit i? ^ we reco ver the vacuum-to- 
vacuum kink (or anti-kink) of Eqn. (^^) with periodic- 
ity /3 ^ 00, corresponding to a zero temperature winding 
transition of the pendulum by quantum tunneling. The 
second limit, E Egph = 2^^ and (3 2tt/u} corre- 
sponds to thermal activation through the unstable, static 
sphaleron solution at q = tt halfway between the two 
neighboring periodic vacuu m sta tes. The action of the 
periodic instanton solution ( 2.10 ) interpolates smoothly 
between these two limits, as illustrated in Fig. 1. 
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Since V"{q = tt) = — the eigenvalues of the sec- 
ond order fluctuation operator of the EucUdean action 
with periodic boundary conditions, expanded around the 
sphaleron solution are 



A. 



27rnY 



(2.14) 



Hence there is clearly exactly one negative mode cor- 
responding to the unstable n = perturbation of the 
sphaleron solution in the high temperature regime (3 < 
/3_, where /3_ is just the period of the harmonic oscilla- 
tion about the inverted potential — V{q) ai q — ir. Ther- 
mal activation through the sphaleron dominates the tran- 
sition rate for high temperatures T > T_ . 

For /3 > /3_, i.e. low temperatures T < T^, the 
sphaleron has more than one negative mode and does 
not dominate the semiclassical winding number transi- 
tion rate. However, in this low temp erature range the 
periodic instanton given by Eqn. ( 2.10 ) possesses exactly 
one negative mode by the following standard argument. 
Since the periodic instanton obeys the classical Euclidean 
Euler-Lagrange equation, 



fq 



V'{q) 



(2.15) 



it follows by differentiation that its time derivative is a 
zero mode of the second order fluctuation operator, i.e. 



(2.16) 



From the periodic behavior of q{T) it follows that q has 
exactly one node in the interval (0, /?), and hence by a 
standard theorem of second order linear differential op- 
erators, the fluctuation operator, 

-TT + = -TT + ^i^^khn^iLUT, k) - uj^ (2.17) 

must possess exactly one mode with no nodes and a lower 
(i.e. negative) eigenvalue. He nce the finite energy pe- 
riodic instanton solution ( 2.10| ) dominates the winding 
number transition rate for temperatures T < Lojlix and 
energies E < 2w^. 

The one negative mode may be understood intuitively 
as the result of the attractive interaction between a kink 
and anti-kink which binds the pair into a kink/anti-kink 
"molecule" and lowers the classical action of the bound 
periodic instanton configuration. This attractive interac- 
tion may be studied analytically in perturbation theory 
at low energies. Since the same approach proves quite 
useful in more complicated models let us review the gen- 
eral method in the pendulum problem. 

Since the kink solution qo{T) of Eqn. (2.6) has zero 



energy we expect the periodic solution for small but finite 
energy to consist of an infinite chain of widely separated 
alternating kinks and anti-kinks which are very loosely 
bound. Thus we consider the trial configuration, 



Qir) 




< T < /3/4 
/?/4 < T < 3/3/4 
3/3/4 <T<(3 



(2.18) 



defined on the fundamental interval [0, /3], and consisting 
of a kink at and /3 and an anti-kink halfway between 
at /3/2. This configuration may be repeated indefinitely 
along the t axis to form an infinite chain of alternating 
kinks and anti-kinks. The configuration is everywhere 
continuous and has finite action per period. 



S[q{T)]= f dT 
Jo 



= 4w^ 



1 

/3/4 



13/4 
3/3/4 



' + Viq) 
dr sech^wT 

f3 



dr sech ui 



/3/4 



= 16cj tanh | — ^ 
V 4 

= 16a; — 32Lue^ 



JI3/2 



0\ 



(2.19) 



where we retain only the leading /3 dependence for large 
p. 



The configurati on ( 2T8) is also a solution of the Eu- 
clidean equation ( ^.15|) except at the matching points, 
T = /3/4 and r — 3/3/4 where the first derivative is dis- 
continuous. That is, 



2(5 T- 



V'iq) 
4 



26 



3/3\ . 

t) 



= Alo sech 



w/3 



-6 T- 



S T- 



3/3 



(2.20) 



These delta functions act as point sources for th e fiel d 
qir). If we expand about the trial configuration (2.18), 
q q + £,, then the non-zero result of Eqn. ( 2.20| ) implies 
that there are linear source terms for the fluctuations ^. 
The effect of these sources on the action functional may 
be found by shifting the fluctuation 
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dT'Gp{r,T')L[q{T')] 



(2.21) 



to remove the linear term ||1 
bmction Gp{t,t') obeys 



|13||. The periodic Green's 



Gp(T,T')^5ti(T~T') 



1 



E 



^27ri(r — r')n//3 



(2.22) 



To leading order in the exponentially small tail e~"^/^ 
we can replace V"{q) by its vacuum value, uP' and use 
the free (perturbative) periodic Green's function. 
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Gp{r,T') 



1 



2cj sinh (/5w/2) 



cosht 



T - r - 



(2.23) 



for r and r' in the fundamental interval [0,/3]. 

To leading order in L[q] the effect of the shift ( ^.21 ) to 
remove the linear term is to alter the action at quadratic 
order, i.e. 



S[q{r) 



S[qir)] 



-\ I dr f dr' L[q[T)]Gp{T,T')L[q{T' 
^ Jo Jo 



)]. (2.24) 



Substituting (2.20) and (2.23) into this extra term leads 



to the result that the action of the periodic instan ton i s 
the action of the trial kink/anti-kink configuration ( 2.1£ ) 
shifted by 



/3 3/3 



r ,'3/3 /3 



r, ,3/3 3/3 



-32t^e-"''/2 + C'(e-"^) . 



(2.25) 



Adding this shift to ( ^.1£ ) we find that the action of the 
periodic instanton in the pendulum model is 



5'(/3) = \&u - 64 we-"'^/^ + O {e-'^'^) 
while its energy becomes 



Em 



dS_ 

dp 



(2.26) 



(2.27) 



to leading order in the expansion in powers of E/E snh = 
16 exp(— tt'/3/2). The negative second term in Eqn. ( 2.26| ) 
is the effect of the short-ranged attractive interaction be- 
tween neighboring kinks and anti-kinks along the chain 
that lowers the action functional to first order in the two- 
body interaction exp(— w(|ri — T2I)) between them. 

This result and this general method of patching to- 
gether zero energy kinks and anti-kinks can be checked 
in the pendulum example by comparison to the exact re- 
sults (p77-2.9) in terms of the elliptic functions. Indeed 
in the low energy limit the modulus A: — > 1 and the ex- 
pansion of the complete elliptic integrals in powers of the 
complementary modulus, 



k' = Vl -fc2 



4e- 



(2.28) 



yields 



5(/3) = 16c.(l- — +o(fc'4; 



= 16 cj — 64 we 



(2.29) 



E{[3) = 2uj^k'^ = 32^2 e-'^/3/2 + q 



(2.30) 



which coincides with (I^H) and ([2!27| ). The exact and 
perturbative results for the action as a function of /3 are 
compared in Fig. 1. 
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FIG. 2. The action of the periodic instanton and sphaleron 
as a function of /3 for the simple pendulum model. The instan- 
ton solution goes from real (solid curve) to complex (dashed 
curve) at /3 = /3_, where the sphaleron action S = EsphP 
(straight line) is tangent to the curve. At /3 = /3cr = 8/lu > P- 
the zero energy instanton action 2So equals the sphaleron ac- 
tion. 

As (3 is decreased from infinity to /3_ and E is increased 
from zero to E^ph, the kink/anti-kink molecule becomes 
more tightly bound and the nonlinear interactions be- 
tween the kink and anti-kink beco me m ore important. 
The action of the periodic instanton ( 2.10 ) decreases from 
25o 16uj to S'(/3_) = EsphP- = 47rtj < 2So and it is 
easy to see from the properties of the elliptic functions 
involved that the derivative, dj3{E) / dE is negative every- 
where in the interval (0, Ssp/i). M E = Egph, /3 = /3_ , 
S = EsphP-, and the periodic instanton solution ( 2.10| ) 
becomes identical to the sphaleron. Beyond this point we 
no longer have a real periodic instanton solution. This 
presents no difficulty as analytic continuation of the solu- 
tion (2.10) for E > Esph is easily accomplished by allow- 
ing the modulus k to become purely imaginary. Then the 
solution q{T) becomes complex (in fact, tt p lus a purely 
imaginary function of r), but the period (2.7) and action 
(p.8[ ) remain real and continue to decrease as the energy 
is increased. However, the fluctuation operator ( 2.17| ) 
now has additional negative modes and hence the analyt- 
ically continued solution no longer dominates the wind- 
ing number transition rate, its place having been taken 
by the static sphaleron for T > w/27r. At /3 = 0.61627/3_ 
and E = 5.74914i?sp;i the action (2.8) of the complex in- 



and 



Stanton actually plunges through zero and then becomes 
negative, finally going to —00 as i? — > 00 and /3 — > 0. 
The situation is most conveniently illustrated graphically 
in Fig. 2. Notice that, in this simple pendulum model 
the zero temperature tunneling exponent 250 is equal to 
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the thermal activation Boltzmann exponent Egph/T at 
T = uj/8 which is less than r_ = uj/2'k. 




plex) direction of the zero mode in function space. Hence 
near /3 = /?_ the periodic instanton solution is given ap- 
proximately by a small perturbation of the sphaleron in 
the relevant zero mode direction, viz. 



FIGS. 3. The action S as a function of the variables 
(a), corresponding to the two near zero mode directions, and 
as a function of one of these variables and the negative mode 
direction (b), for fixed /3 slightly larger than /3_. The 
ring of extrema of S in the first figure (due to the arbitrary 
phase in Eqn. 2.3f ) or the two outer extrema in the second 



figure correspond to the periodic instanton solution. As /? ap- 
proaches (3- these extrema merge with the sphaleron in the 
center. Since the periodic instanton solution has one nega- 
tive mode and lower action than the sphaleron, it dominates 
the winding number transition rate which is controlled by the 
lowest available saddle in (b). 

The merging and resplitting off into the complex do- 
main of the periodic instanton solution as (3 is decreased 
through /?_ is an interesting aspect of the behavior we 
have just sketched. This may be understood from the 
fact that near (3 = [3- the eigenvalues of the static 
sphaleron fluctuation operator, A±i in ( p. 14 ), go through 
zero. Whenever a zero mode appears in the second or- 
der fluctuation operator this indicates the existence of a 
nearby solution of the classical equations. When the zero 
mode is not related to a symmetry of the action but in- 
stead appears only at special values of a parameter in the 
boundary conditions, then the solution generally splits off 
or bifurcates into two (or more) different solutions as the 
parameter is varied. Conversely, as the parameter is var- 
ied in the opposite direction two or more solutions merge 
at the critical value of the parameter. This is just the 
case for the static sphaleron solution as the parameter 
(3 is varied through /9_. At (3- we find another classical 
solution with the same Euclidean period splitting off or 
merging with the static sphaleron in the (generally com- 



q{T) ~ Qsph + 6 cos{lut + . 



(2.31) 



where (fi is an arbitary phase of the periodic solution and 
5 is a small parameter that goes to zero as /? — *■ /?_. 
The arbitary phase informs us that in fact there is a 
one parameter U{1) family of periodic instanton solutions 
which merge and then resplit off from the sphaleron at /3_ 
which accounts for the two zero mode directions, n = ±1 
at /3 = 

This situation is also easier to visualize graphically. Let 
us suppress all but the relevant zero mode directions and 
represent the action functional 5" as a function only of a 
few variables £,i together with the parameter (3. Clearly 
as f3 varies the surfaces oi S = const, change, as do the 
location of the extrema, dS/d£^i = 0. At certain values 
of P if two of the extrema split off or merge, there is a 
zero eigenvector of the second order fluctuation matrix 

S / d^id£,j . Since the number of extrema is determined 
by the global properties of S, generally the bifurcation of 
extrema does not mean that their total number changes 
discontinuously, but rather that real extrema move into 
the complex domain as the parameter (3 is varied through 
/3_ . The U{1) invariance of the action under the arbitrary 
phase (f) implies that a ring of equal action extrema con- 
verges on the sphaleron and then becomes complex for 
(3 < (3-. This action as a function of the three relevant 
variables and is illustrated in Figs. 3. 

To summarize this discussion, the physical interpreta- 
tion of the properties of the classical Euclidean periodic 
solutions for winding number transitions in the simple 
pendulum model is clear. At low temperatures T < l//3_ 
the periodic instanton solution dominates the vacuum-to- 
vacuum winding number transition amplitude which is 
exponentially suppressed. As the temperature is raised 
the classical periodic solutions become a more closely 
spaced chain of kinks and anti-kinks with a larger non- 
linear interaction between them that lowers the action. 
The tunneling rate becomes less suppressed as the turn- 
ing points of the periodic solution move inwards toward 
their midpoint aX q — -k, corresponding to less quantum 
tunneling and more thermal motion in the allowed region. 
At T = T_ = ujl-K the turning points become identi- 
cal with the sphaleron at g = tt and tunneling has been 
replaced by purely classical thermal activation which is 
still only slightly less exponentially suppressed than the 
zero temperature transition rate. At this and all higher 
temperatures the winding number transition rate is dom- 
inated by the sphaleron which has exactly one negative 
mode for T > T_ . The periodic instanton solution moves 
off into the complex domain and no longer contributes 
to the transition rate, which becomes less and less expo- 
nentially suppressed as the temperature is raised further. 
Finally at temperatures greater than of order Es^h = 
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the exponential suppression disappears entirely and the 
pendulum swings freely around its pivot, the probability 
of finding the pendulum anywhere then becoming nearly 
uniform as the temperature is raised still further. 

This fixed temperature discussion has a natural analog 
at fixed energy. The probability of making a winding 
number transition at fixed energy between E and E + dE 
may be expressed in the form, 



PiE)dE^J2\(f\^^E\'^\' 



(2.32) 



where S is the 5-niatrix, Ve is a projector onto energy E 
and the initial and final states |i) and |/) lie in different 
winding sectors in the energy interval E to E + dE. Peri- 
odic instantons appear again as the configurations which 
saturate P{E), and the probability is given, with expo- 
nential accuracy, by the exact analog of the quantum- 
mechanical formula, 



where 



P{E)dE - e^]i{-W{E))dE , 



W{E) = S{(3{E)) ~ E(3(E) 



(2.33) 



(2.34) 



is the Legendre transform of S{(3). The relevant instan- 
ton solutions should be considered now as fixed in energy 
rather than in periodicity /3. 

The initial and final multi-particle states can be read 
off from the analytic continuation of the finite energy 
periodic instanton into Minkowski time at its turning 
points. Hence periodic instanton solutions to the clas- 
sical Euclidean equations contain non-trivial information 
about multi-particle transition amplitudes between dif- 
ferent winding number sectors at finite energy, ft has 
been suggested that by suitably modifying the boundary 
conditions in the complex time plane, information about 
transition amplitudes from few particle initial states to 
many particle final states may be obtained as well |ll0| . 

The fixed temperature transition rate can be recon- 
structed from the fixed energy tunneling probability by 
weighting with a Boltzmann distribution and integrating 
over energy. To exponential accuracy we have 

r(r)-^°°di?exp(^-|^p(i?) 

dE exp(^-^-W{E)^ . (2.35) 

Evaluating this integral by the method of steepest de- 
scent (which is valid in the limit of arbitrarily weak cou- 
pling uj~^ ~ — > 0, provided T is coupling independent) 
gives the saddle point condition. 



1 

T 



dW 
~dE 



I3{E). 



(2.36) 



by the properties of the Legendre transform W = S—E(3. 
Hence we recover the connection between the tempera- 
ture and Euclidean periodicity of the instanton solution 



in this wa y. Ev aluating the second derivative of the ex- 
ponent in ( 2.35 ) at this saddle point yields 



' dE"^ 



dP{E) 
dE 



(2.37) 



wh ich te lls us that only solutions for which the quantity 
in (2.37) is negative can contribute to the fixed tempera- 
ture rate r(T). On the other hand, classical solutions for 
which thi s der ivative is positive cannot contribute to the 
integral (2.35) since they are local minima rather than 
maxima of the the exponent in the integrand. In the case 
of the simple quantum pendulum example considered in 
detail in this section ( 2.37] ) i s ind eed negative and the 
periodic instantons given by ( 2.f 0| ) do contribute to the 
finite temperature transition rate T{T) in the expected 
way, for all temperatures T < r_. We turn now to the 
0(3) sigma model where the behavior of the p eriod ic in- 
stanton solutions, the corresponding quantity (2.37) and 
the crossover from quantum tunneling to thermal activa- 
tion are all quite different. 



III. THE 0(3) SIGMA MODEL 

The two-dimensional 0(3) nonlinear sigma model is 
defined by the Euclidean action functional. 



5* - ^ / dxdrid^ny , 



(3.1) 



where n'^{x), a — 1,2,3 are three components of a unit 
vector, n°"n°' — 1, and fi — t,x. The constraint on the 
magnitude of n° at every spacetime point is the source 
of nonlinearity in the model. A convenient parameteri- 
zation in which the constraint is removed at the price of 
explicit nonlinearity is the complex field definition. 



ni + in2 
w = . 

1-713 

Introducing the complex coordinate 
z = X + IT = re 



(3.2) 



(3.3) 



in terms of the spatial position x and Euclidean time t 
enables us to rewrite the action ( |3.l|) in the equivalent 
forms, 



S=^[dxdr^^^ 
J (1 + ww)'^ 



S = 



9^ 



dx dr ■ 



1 



or 

dw dw dw dw 



(1 -I- ww)'^ \ dz dz dz dz 



(3.4) 



where the overbar denotes complex conjugation and w — 
w{z, z) — w{z, z). 

This model possesses some remarkable similarities with 
pure gauge theories in four dimensions. The most impor- 
tant properties which concern us here are the conformal 
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invariance of the classical action S, and the existence of 
a periodic vacuum structure in the model. In connection 
with the conformal invariance we observe that the cou- 
pling constant g is dimensionless and there is no length 
scale in S. Correspondingly, the quantum theory is log- 
arithmically renormalizable and in fact, asymptotically 
free in the coupling g . 

The periodic structure should be clear from the fact 
that an 0(3) rotation by 2tt around any axis brings the 
vector n° back to itself. The topological winding number 
associated with this field periodicity will be made explicit 
if we identify the points at infinity of the complex plane 
in the coordinate z — x + it. Then the plane has the 
topology of the sphere S^. Since the field n'^ is also con- 
strained to lie on 5^, the n° field is a map from 5^ to 
5^ which is characterized by an integer winding number, 
given explicitly by 



1 

8^ 
1 

TT 



dx dr 



1 / dw dw dw dw 
(1 + wwY \ dz dz dz dz 



(3.5) 



Comparing the latter form of the topological winding 
number with the action (3.4) it is clear that the action in 



any integer Q topological sector is bounded from below, 
i.e. 



S> 



An 

7 



(3.6) 



and moreover that this bound is saturated by meromor- 
phic functions w of the complex variable z or z ]l5| . In 
particular the conformal map of to 5^, topologically 
equivalent to the identity map is the one instanton solu- 
tion of the Euclidean equations with Q = 1. This solution 
can be written in the form, 



wo{z) = - = -e 
z r 



■ie 



and has action 



47r 

•^o — ^ 
9^ 



(3.7) 



(3.8) 



The anti-instanton solution is obtained from (3.7) by Eu- 



clidean time reversal which in the w description is equiv- 
alent to complex conjugation. 

The instanton and anti-instanton solution have zero 
Euclidean energy. 



E^^^ ( dx i^drn^drn" + d^n^d^n") 
^9 J 



0. (3.9) 



They correspond respectively to the classical Euclidean 
path of least action and its time reverse, connecting peri- 
odic vacua separated by unit winding number. Together 
the instanton and anti-instanton constitute a periodic 
trajectory beginning and ending at the same vacuum 



state, with one negative mode corresponding to the at- 
tractive interaction between them, just as in the simple 
pendulum model. Hence in the limit of weak coupling 
g^ — > the rate for a unit winding number vacuum to vac- 
uum tunneling transition is exp(— 2S'o) = exp(— 87r/5^) to 
exponential accuracy, and is strongly suppressed. 

Because of the conformal invariance of the action (^.l|), 
the zero energy instanton can have arbitrary scale p. This 
implies that the potential energy barrier between winding 
number sectors, i.e. the second term of ( |3.9| ) evaluated 
at r = which is proportional to can be made 

arbitrarily small and no finite energy sphaleron solution 
exists in the symmetric 0(3) model. This feature is quite 
different from the simple pendulum and abelian Higgs 
models but is shared by the pure Yang-Mills theory in 
four dimensions. 

In the electroweak theory conformal invariance is bro- 
ken in the Higgs sector which then makes possible the 
existence of a classical sphaleron solution. In the 0(3) 
sigma model the conf orm al invariance may be broken by 
adding to the action (B.l) the explicit mass term 



(3.10) 



which also violates the 0(3) symmetry and fixes the vac- 
uum state to be n^^yg^^-^ = (0,0,-1). In the Heisenberg 
spin language this corresponds to placing the system in 
an external magnetic field which aligns all the spins in 
the direction n'^^^^-^ at zero temperature. 

When the sigma model has been modified in this way 
by the introduction of the soft symmetry breaking mass 
term 5^, 



s + s 

77 



(3.11) 



the instanton configuration with arbitary non-zero p in 
( |3.7|) ceases to be an exact solution of the field equations, 
but there is now an exact sphaleron solution | p^ , p7t , 
namely. 



(sph) 
fsph)ix) = 
fsph) (^) = ~ 



2 tanh(mx) sech(TOa;) , 

1 -h 2sech^(mx) (3.12) 



with finite energy, 



^sph 



8m 



(3.13) 



Geometrically, this solution maps the infinite spatial line 
onto a great circle beginning and ending at the south 
pole, fi^^^^-j = —1. The sphaleron lies exactly halfway 
between the vacua of winding number and 1, and Esph 
is therefore the height of the potential energy barrier be- 
tween these vacuum states, which is now fixed and finite. 

Let us recall that the spectrum of time independent 
perturbations around the static sphaleron has exactly 



8 



one normalizable negative mode with the eigenfunction 
ut{x) = (0, sech^(mx), 0) and the eigenvalue e?_ = 
— 3m^. There are two zero modes of the spatial fluctua- 
tion operator corresponding to spatial translations of the 
sphaleron position and rotation of its great circle trajec- 
tory around the axis. The remaining eigenvalues are 
strictly positive and form a continuous spectrum with 
e+ > . Taking the periodic r dependence of the fluc- 
tuations into account yields the following eigenvalues for 
the full second order fluctuation operator of the action 
around the sphaleron solution, 



An 



(3.14) 



Periodic boundary conditions are the appropriate ones 
for contributions to the path integral at finite temper- 
ature T = p-'^. By the same argument as in the one 
dimensional pendulu m th e existence of one and only one 
negative eigenvalue ( 3.14 ), namely n = and = — 
-3m2 for (3 < (3^ = 27r/|e_| = 2n/{V3m) = 3.628 m"^ 
implies that the sphaleron configuration contributes to 
the winding number transition rate for temperatures 
T > T_. The sphaleron contribution to the rate per 
unit one dimensional spatial volume can be expressed in 
the notation of the present paper in the form, n& 



^ sph 



where 

h{m/T) 
4m 
vr 



mT 



7rg2 sin(|e_|/2r) 



e:>q)[-Esph/T - 



h{m/T)] , 

(3.15) 



dk 







-I- 



In (l - e-'^'-l'^^ , (3.16) 



and 



UJk 



(3.17) 



The asymptotic form of this function in the high temper- 
ature limit m/3 ^ 1 is 

4m 

him/T) ^ -3 ln(m/r) ln(m/T) - C + 0{mlT) , 

ttT 

(3.18) 

where we take this opportunity to correct an error in the 
numerical evaluation of the constant C in Eqn. (5.7) of 
Ref. PI, viz., pi 



C 



dx 



1 



1 



+ 1 x^ +4 



lii{x^ + 1) 



In 12 = 2.4849 



(3.19) 



As a result, the constant defined in Eqn. (5.17) of Ref. 
p6[ changes as well, becoming K = 2e'^/^/3 = 8V3 ~ 
13.8564 . 



The sphaleron thermal activation rate (3.15) is less 
than the instanton mediated tunneling rate unless T > 
Tcr where the crossover temperature. 



T„ = ->T^ = ^ 



(3.20) 



is defined by equating the quantum tunneling and ther- 
mal activation exponents, 25*0 — Esph/Tcr- Hence al- 
though the sphaleron provides a second semiclassical es- 
cape path for winding number transitions for all temper- 
atures above T_, its contribution is exponentially sub- 
dominant until the temperature T > Tcr > where 
there is a sharp crossover from quantum tunneling me- 
diated by the zero scale singular instanton to thermal 
activation via the static sphaleron, in the semiclassical 
weak coupling approximation <^1. 

The fact that the crossover between winding number 
transitions dominated by quantum tunneling and those 
dominated by thermal activation takes place at a tem- 
perature T cr gr eater than that at which the n = 1 



eigenvalue ( 3.14 ) goes through zero is a qualitative dif- 
ference from the situation in one dimensional quantum 
mechanics or the abelian Higgs model in 1 + 1 dimen- 
sions. In each of these examples the temperature at 
which the sphaleron and instanton/anti-instanton expo- 
nents are equal is smaller than the corresponding T_ of 
the sphaleron negative mode. In such cases the crossover 
between zero temperature quantum tunneling and fi- 
nite temperature thermal activation is smooth and the 
sphaleron dominates the transition rate for all T > T_. 
The situation for the 0(3) sigma model is illustrated in 
Fig. 4 where the numerical results from Section 5 are 
incorporated. 




FIG. 4. The action of the periodic instanton and sphaleron 
as a function of P for the 0(3) sigma model. As in Fig. 2 
the instanton solution goes from real (solid curve) to complex 
(dot-dashed curve) at /3 = /?_, where the sphaleron action 
S — EsphP is tangent to the curve. However, in contrast to 
Fig. 2 the value of Per at which the sphaleron action equals 
the zero energy action 2^0 (dashed horizontal line) is less than 

As is easy to see from scaling arguments (Derrick's 
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theorem) the mass term 5*^ now drives the instanton size 
p — > 0, so there is no exact instanton of the form (3.7) 
except the singular one with p ^ 0. Ho wev er, for small 
but finite p the instanton configuration (3/7) remains an 
approximate solution for r ^ m^^. Furthermore, for 
r ^ p, I Wo I ^ 1 and the field equations linearize for any 
TO. Hence if mp <C 1 a smooth field configuration which 
is an approximate solution for all r is 



i{x,T)=mpe ^ Ki{mr) , 



(3.21) 



where the Bessel function Ki is the solution of the lin- 
earized Euler-Lagrange equations obeying the boundary 
condition at infinity. This approximate solution with ac- 
tion close to the instanton action An/g^ contributes to 
the winding number vacuum to va cuum transition am- 
plitude just as the exact instanton (3^) solution does in 
the symmetric m = model. 

In either the softly broken or unbroken sigma model 
instanton/anti-instanton configurations are also approx- 
imate solutions to the field equations with finite energy. 
In the low energy regime perturbation theory is valid and 
one can construct the approximate solutions by form- 
ing an infinite chain of alternating instantons and anti- 
instantons along the Euclidean time axis separated by a 
half period /3/2. At this point a new parameter, the pe- 
riod (3 enters the problem and we have a finite interaction 
between instantons and anti-instantons along the chain, 
just as in the pendulum example of the previous section. 
Hence we can have a competition between the tendency 
of the scale p to shrink to zero for an instanton in iso- 
lation and the tendency of p to increase due to the at- 
tractive interaction between neighboring instantons and 
anti-instantons, with the two effects balanced at a par- 
ticular p which is a function of [3. In order to obtain 
some detailed understanding of this qualitative picture 
and provide a quantitative benchmark for the numerical 
methods to follow, we work out in the next section the 
periodic instanton solution for energies E <C E^ph where 
perturbative methods are applicable. 



IV. LOW ENERGY PERTURBATION THEORY 

Our construction of the periodic instanton solution at 
low energies will take place in two steps. Since we are 
not in possession of an exact solution to the classical 
equations f or fin ite mp we evaluate the action of the con- 
figuration (3.21) to first non- vanishing order in mp, at 
first without regard to the finite periodicity (3. Then 
we proceed to patch the configurations together to con- 
struct a trial configuration periodic in Euclidean time as 
in the pendulum example. Extremizing the resulting ac- 
tion function with respect to the free scale parameter p 
will determine the particular p{(3) at which the attractive 
interaction between instantons and anti-instantons just 
balances the self-interaction which would cause an in- 
stanton to shrink to zero size in isolation. Provided that 



in the end mp{f3) <C 1, the expansion is consistent and 
we may trust the resulting extremal action constructed 
in this way as the correct periodic instanton action to 
first non- vanishing order in the small parameter mp{j3), 
or cquivalently E/Esph.. 

The sigma model action for any function of the form 



w{x, t) 



-ie 



f{r) may be expressed as 
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■ dr 



{1 + Pf il + P) 



(4.1) 



For the trial configuration / = mpKi{mr) we make use 
of the Bessel function identity, K[ + Ki/u = —Kg to 
rewrite the first integral in the preceding expression in 
the form. 



2 2 
TO p 



u du 



{Kl{u)-2KiK'Ju) 
[l + m^p^Kl{u)f 



(4.2) 



Let us divide the integration over u = mr into two 
parts, u < a and u > a where a is chosen so that mp <C 
a <C 1 in the limit of small mp. For u < a we can make 
use of the series expansions of the Bessel functions for 
small argument. 



Ko{u) 






Ki{u) 


1 

u 


u 
^2 


K2{u) 


2 


1 

^ 2 



In 



+ IE 
''0('„,2 1 



IE 



(4.3) 



taking care to include all terms that give a contribution 
to the final answer up to order m^ p^ , while in the interval 
u> a we can safely replace the denominator in (4.2) by 
unity since the terms neglected are higher order in m^p^ 
and the integral with lower cutoff a is non-singular. In 
this way we find 



h 



u du 



[l + m^p^Kf{u)f 



7B + o 



(4.4) 



+0{m p In mp, a In a) 



and 

l2 = 



I du 



{Kl{u)-2KiK'Ju) 



[l + m^p^Kl{u)Y 
= ^ + In +1E + 0{m^p'^ \n m.p, \n a) (4.5) 

where — 0.577 ... is Eulcr's constant. In the lat- 
ter integral we have made use of the integration formula 
5.54(2) of Ref. pl. 



u duKQ{u) 



a 

y 



[K',ia)-K',{a)] . (4.6) 
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Combining the results (4.4) and (|4.5|), the a dependence 
drops out (as it must) and we obtain 



h+l2 = 



1 



0{m p In mp) 



(4.7) 



Analyzing the second integral in (4A) in the same way 
we obtain 



udu Ki{u) 



and 

h = 



[l + mV^?(")] 
= Ina — In 771/9 + 0{m^ p^ In^ mp, h\ a) (4.8) 



udu K1{u) 



[l + m^p^K^{u)] 
= - In (^) - 7B - i + 0{m'^p^ In^ mp, \i? a) (4.9) 

where in this latter integral we have made use of the 
integration formula, 



uduKliu) = — [Ko{a)K2{a) - Kl{a)\ . (4.10) 



Hence we find 

h + h = - 



In {-if) - - ^ + 0{m^p^ \n mp) , 



(4.11) 



and all together the action of the trial configuration is 
47r?7i^p^ 



47r 
9 



(4.12) 



In (^) + 7b] + 0(mV' In' mp)} 



The leading 'in/g^ term is clearly the one instanton ac- 
tion of the unbroken [m = 0) model, so that in con- 
sidering periodic instantons next we should multiply the 
result ( 4.12| ) by two. 

In order to calculate the corrections to the action due 
to the finite periodicity in the sigma model we follow the 
same method as in the pendulum example. We consider 
the periodic configuration. 



u>(r, a:) 



Wm(T,x) , 



< r < /3/4 
/3/4 < T < 3/3/4 
3/3/4 <T<(3 



(4.13) 



and evaluate its action in the fundamental interval [0, /3] 
to first non-trivial order in m^p^. Starting from 



r Jo 



P roo 

dT 



(l + |u;|2)2 



w ■ 



(4.14) 



the /3 dependence of the action comes from the tail of 
the configuration where \w\ <C 1, and hence to leading 
order in m^p^ the jwp terms in the denominators may be 
neglected. By introducing the spatial Fourier transform 



dx 



Wm{T, k) 

—2imp 
-i — {ujk ~ A:)exp(-|T|wi;) 



dxe WmiT,x) 



(t cos kx — X sin kx) Ki (mr) 



(4.15) 



we find that the /3 dependent part of S may be expressed 
in the form, 

dS 8 d fP'\ r dk r,^ ^ . ,.|2 
df3 52 271 V I 

^ir,k)\'} 



47rp2 



9 Ja 



dk {ujI + fc2) exp 



(4.16) 



to leading order in m^p^. The infinite periodicity, /3 
jTidependent part of S cannot be evaluated from only 
the quadratic terms in the action but it is just given 
by twice the zero energy action we calculated in (4.12) 
above. Hence, 



5N-^{l-mV 



ln(^)+7^. 



r,fcM±^exp 



' ^2^ "^ + 0(m'*p* In' mp) 

(4.17) 



is the action of the configuration ( 4.13 ) to this order. 

As in the pendulum example of the last section the 
joining of the configurations at r = /3/4 and t = 3/3/4 
produces delta function source terms proportional to 

r, ^ f 9^ 9^ 2\ , ^ 

^\^o'^ ( ^ 3/3 
-2drWrn(-,x)d I ^ ^ I + 2dr (- - , a;)d I - — 

(4.18) 

Making use of the periodic thermal propagator in two 
dimensions. 



Gp{t,t';x,x') = 
dk 1 



_^ 2tt 2u;k sinh {(3uJk/'2) 



coshtjfe \t — t' 



(4.19) 



and taking account of the two real degrees of freedom in 
the complex field w, we obtain the quadratic shift in the 
action, 
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^ J dr J dr' J dx J dx'L'[w*]Gp{T, t' ; x, x')L[w] . 

(4.20) 



Substituting ( 4.18| ) and ( 4.19 ) into this expression, we 
find that the quadratic shift becomes 



dk 

T^^k 
1 



\drWr, 



— , fc I 1^ coth ( 

4' ; I \ 2 ) 



sinh(;3c.,/2T^^<U'^'^^"" 



(4.21) 



We did not need to keep the analog of the second term 
above in the pendulum model since there it is subdomi- 
nant to the first term as /3 — > oo. However, in the 0(3) 
model /3 ^ as ^ as we shall see presently, and this 
second term must b e retained as well. 
Since from ( 4.15| ), 



d^w„ ( A: ) = iTrp (wfc =F k) exp 



(4.22) 



the q uadra tic correction to the action arising from the 
shift ( |420| ) is 



dk 



1 



g2 Wfc sinh(/3wfe/2) 

^2 



k' 1 



) 



(1 



(4.23) 



Combining this shift with the unshifted action (4.17) we 
secure finally. 



5(/3,p) = ^{l 



2 2 

m p 



i^(m/3)+ In (^)+7B 



to the lowest non-trivial order in mp, where 



F(mP) 



dk 



1 



Wfc sinh(/3a;fe/2) 



2fc2 



(4.24) 



> 0. 



(4.25) 



This result may also be obtained by the "R-term" method 

With the action of the trial configuration in hand for 
arbitrary small mp we can determine the value of p which 
extremizes the action and therefore leads to a classical 
periodic instanton solution close to the trial configuration 
(p3|). Since 



dS{p,p) 
dp 



vanishes for p = or for 



2F + 21n(^)+27i, + l 



(4.26) 



P = PiP) 



exp -'F{m(3) - 7_e - 



(4.27) 



the zero energy singular instanton at p = with finite 
action S = 877/(7^ is always a solution of the equations 
for any period (3, but there is also a non-trivial periodic 
instanton solution with action, 



5(/3) = 5(/3,p(/3)) = ^ 
9^ 



and energy. 



1 + 2e"^^("'')"^'^ 



E-l 



> 



Stt 



9^ 
(4.28) 



327rTO 



g-2F(m/3)-27.-lj^/(^^) 



-2F 



167re-2T^-i 

2 ' 

9 

— ^ -I- 1 I cosh 



dk 



sinh2(cjfe/3/2) 
2 



1 



(4.29) 



This second non-singular solution with non-zero p(/3) is 
the periodic instanton solution we seek in the limit of 
small but finite E. 

The qualitative picture sketched in the previous section 
has been verified by this explicit calcu lation, namely at 
the non-trivial value of p(/3) in ( 4.27 ) the tendency of 
an individual instanton or anti-instant on to contract to 
p = 0, represented by the logarithm in (4.26) is just bal- 



anced by the attractive interaction between members of 
the infinite chain which tends to increase p, represented 
by the F{m/3) term in (4.26). This competition between 
the two effects can occur in the sigma model only be- 
cause of the softly broken conformal invariance of the 
mass term 5*^ which leads to the existence of the confor- 
mal scale parameter p. 

Notice also that in contrast to the pendulum example 
of the last section it is not /3 — > 00 which characterizes 
the low energy limit E 0, but rather F{m(3) — > cxd in 
the sigma model. Since F{m(3) vanishes exponentially as 
/3 — > 00 but in the opposite limit. 



F(m(3) 



Tin?fp 



00 



m/3 ^ , (4.30) 



the low energy limit is characterized hy P 0. This 
(possibly counterintuitive) result is nevertheless consis- 
tent with perturbation theory and the dilute gas picture 
of instantons and anti-instantons at low energies because 
the size of the individual (anti-)instantons goes to zero 
much more rapidly with /3 — > 0, and it is the ratio, 



m 2£ 

p 



■ exp 







(4.31) 



which controls the validity of the dilute gas approxima- 
tion. Clearly then as /9 0, p(/3) — > exponentially, 
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an d it is legitimate to neglect the higher powers of mp 
in ( 4.29 ) in the low energy limit, E <C Egpin which jus- 
tifies our expansion in powers of mp a posteriori. In the 
extreme low energy limit we have the asymptotic forms. 



S{(i) 



647r3 



g2g27_B + lTO2/33 



■ exp 



exp 



27r" 
rrfip'^ 



and 
(4.32) 



From either this last expression or the more accurate 
(4.29) above we find that the behavior of E increasing 
with increasing (3, 



dE 
'dp 



> 



(4.33) 



is actually the case for the periodic instantons in the low 
energy limit of the nonlinear sigma model. We shall see 
in the next section that this behavior persists for all ener- 
gies up to the sphaleron energy. Hence it is already clear 
that the behavior of the periodic instantons in the 0(3) 
model is very different from that of systems with only 
one quantum mechanical degree of freedom such as the 
simple pendulum, or the abelian Higgs model. In those 
models f3{E) ^ oo as -E ^ 0, and the periodic instan- 
ton contributes to the transition rate as the temperature 
goes smoothly to zero. From the discussion at the 
end of the second section it follows that these periodic 
instantons in the sigma model with period going to zero 
as -E — > cannot contribute to the low temperature tran- 
sition rate, although they can contribute to, and in fact 
dominate, transitions between noK-thermal low energy 
states. 

Clearly this behavior is a direct consequence of the con- 
formal invariance of the unbroken model which leads to 
the additional scale parameter p in the instanton config- 
uration, and as it turns out, an additional negative mode 
direction in the fluctuations around this conflguration. 
This additional negative mode is easily checked from the 
second derivative of the classical action (4.24), 



dp^ 



167r 



(3 + 2jE 



-2 In 



2F(m(3) 



-TO' 



< 



P=p(/9) 

(4.34) 



on the solution with non-zero p(/3). Thus, any small 
change of p from its critical value, p{(3) decreases the 
action of the configuration and implies the existence of a 
negative eigenvalue in the second order fluctuation oper- 
ator around the non-singular periodic instanton solution 
that is quite independent of the usual negative mode with 
fixed p as found by the standard argument of the Section 
2. The origin of this second negative mode about the 
non-trivial periodic instanton solution we have found for 



low E should be clear from what has already been said 
regarding the balance between the attractive interaction 
between nearest neighbor instantons and anti-instantons 
on the one hand, and the self-interaction attracting each 
individual (anti-)instanton towards zero size on the other. 
If the precise balance between the two is upset by vary- 
ing p slightly away from its critical value, /o(/3), then one 
or the other of the two interactions dominates and p is 
driven further away from its critical value, i. e. there is a 
negative mode in the direction of varying p, which is just 
what (4.34) makes explicit. 




FIGS. 5. The action S plotted as a function of the same 
variables ^±1 (a), and (b), for fixed /3 < as in Figs. 
3, but for the 0(3) sigma model. The periodic instanton is 
again represented by the ring of extrema of S in (a) or the two 
outer extrema in (b), which merge with the sphaleron in the 
center as /3 approaches /9_ . Since the periodic instanton solu- 
tion now has two negative modes and higher action than the 
sphaleron, it is subdominant. The finite temperature wind- 
ing number transition goes over the lowest saddle, either the 
sphaleron or the p = instanton (not shown) depending on 
the temperature. 

It is this second negative mode direction due to confor- 
mal rescalings that makes the non-trivial periodic instan- 
ton configuration subdominant at finite temperature, at 
least in the low energy limit where perturbation theory 
and the form ( 4.29| ) a pplies . Correspondingly, the peri- 
odic instanton action ( 4.28| ) is greater than the singular 
solution with p = 0. Near the sphaleron energy the ac- 
tion as a function of the three relevant variables ^±1 and 
^0 for the sigma model is illustrat ed in Figs. 5. 

From the first form of Eqn. (4.34) we observe that 
the zero size instanton/anti-instanton configuration has 
no second negative mode in the p direction. It has only 



13 



the expected single negative mode corresponding to the 
attractive interaction between the pair. Hence this sin- 
gular configuration with fixed finite action per period 
2So (shown as the constant straight line in Fig. 4) can 
and does contribute to the winding number transition 
rate at low temperature T < Tcr- Notice from Fig. 4 
that although the actions 25o and EsphP are equal at 
(3=1 /Tcr the actual field configurations of the zero size 
instanton/anti-instanton pair and the sphaleron are quite 
different and there is no merging of these solutions at 

The picture we have sketched in this section is based on 
the existence of a perturbative instanton/anti-instanton 
expansion at low energies and knowledge of the sphaleron 
and its negative mode at finite energy and (3 near [3cr- 
Continuity of the space of solutions to the classical Eu- 
clidean equations suggests that the periodic instanton so- 
lutions merge with the sphaleron at [3 approaches Per 
from below with dE/df3 > 0. However one does not 
know how the action and energy of these solutions will 
vary as P is varied from f3cr « priori. In the absence of 
analytic methods for finding the periodic instanton so- 
lutions to the Euclidean equations in the intermediate 
region between _B — > where perturbation theory holds 
and E Esph where the solutions approach the static 
sphaleron, one must rely on a numerical approach. It is 
to the details and results of this numerical study which 
justifies Figs. 4 and 5 that we turn next. 

V. NUMERICAL TECHNIQUE AND RESULTS 



then returns to its initial configuration at r = /3 by sim- 
ply reversing the sign of all r derivatives. Thus we enforce 
the half-period boundary conditions, 

dn-{T = 0,x) dn-{T = f3/2,x) 

d-r = d-r = '^' ^^-^^ 

which removes the time translation invariance of the so- 
lution. The other boundary condition we require is that 
at spatial infinity the solution approaches the vacuum. 
Since we shall work in a finite box of length L we require 

n'^ir, X = -L/2) = n"(r, x - L/2) = = (0, 0, -1) 

(5.5) 

for mL large but finite. 







X 
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FIG. 6. The finite rectangular region (r, a;) of coordinate 
space that maps onto the shaded region of the sphere in Fig. 
7. 



In this section we present the results of a numerical 
study of periodic instanton solutions in the two dimen- 
sional nonlinear 0(3) sigma model modified by the mass 
term Sra ■ A preliminary version of these results has been 
reported earlier in Rcf. At energies comparable, but 
not very close to the sphaleron energy, the periodic in- 
stanton solution has to be found numerically, i.e. one 
must solve the Euclidean field equations, 

-{dl + dl)n'' ^m^bl-an" = (5.1) 

together with the constraint, 

n^n" = 1 . (5.2) 

Here, a is a Lagrange multiplier enfo rcin g the constraint, 
whic h is easily found by multiplying ( ^.l[ ) by n° and using 
(^.2|), namely 

= df.n^df.n" + m^ns . (5.3) 

The periodic instanton solution we seek has vanishing 
time derivative at initial time t — 0, evolves to another 
turning point at half-period /3/2 where it reflects and 




FIG. 7. Geometric representation of the periodic instan- 
ton range on S^. 

Geometrically, the periodic instanton solution obeying 
these boundary conditions maps the rectangular region 
of (r, x) coordinate space pictured in Fig. 6 into the 
shaded region of the sphere in Fig. 7. From these fig- 
ures it should be clear that the solution may be chosen 
to have well-defined symmetry properties under reflec- 
tion through the lines bisecting the rectangular region of 
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Fig. 6, so that at a; = 0, rii vanishes while n2 and 713 
reach extrema, and at r = /3/4, n2 vanishes while ni and 
n3 reach extrema. Imposing these symmetries fixes com- 
pletely the spatial translational invariance and rotation 
invariance around the 713 axis of the solution, as well as 
reduces the region we must consider to only one quarter 
of the full rec tan gle in Fig. 6. Hence we solve the differ- 
ential Eqns. (5.1) subject to the boundary conditions, 



X = -L/2 
T = 

T = f3/4 



rii =0; = 0; nj, -1; 
drill dr = dn2/dT = dn^/dr = 0; 
ni = 0; dn2/dx — dn^/dx = 0; 
n2 — 0; dni/dr — dn^/dr — 0. 



(5.6) 



The lattice version of Eqns. {5A) can be obtained by a 
variational principle starting from the discretized action. 
Let {Ti,Xj), where z = 0, . . . , / and j = 0, . . . , J, be the 
coordin ates of lattice sites. The lattice version of the 
action ( ^.l[ ) and ( 3.10 ) reads 



r 1 



dxj 



1 



i,3 



dn 



dxj 



-f 777^ (1 -I- 77^^ )dfidXj I 

where we have used the notations dri — r^+i - 

(l/2)(dTi_i + dTi) and similarly for dxj, dxj. 

The equations which follow from the action (5/7) can 
be written in the form 



(5.7) 



dfi 



[S'^-77'^(B^77b)]y -0, 



(5.8) 



where 



dTi dTi dTidTi. 



dxjdxj-i 



(5.9) 



Eqn. (5.S) actually contains only two independent equa- 



tions since its projection onto the vector nfj is iden tically 



zero. Together with the constraint equation ( |5.2| ) these 
comprise a complete set of three independent e qua tions 
associated with each lattice site. From Eqns. ( ^.8[ ) and 
( 5_^) it is clear that i?" must be anti-parallel to 71" at each 
point (i, .1 )- Since 71" is normalized to unity, the Eqns. 
(5.8)~( |5.9| ) can be rewritten in the equivalent, symmetric 
form. 



(5.10) 



for every (i, j), and where the negative sign of the square 
root has to be taken because UaB"^ < 0. This form is 
most convenient for numerical calcul ations . 

To obtain a numerical solution of ( |3.10| ) we use New- 
ton's method |2l]. That is, we 

1. choose an initial field configuration as a first guess; 



2. linearize the equations (5.10) in the background of 
the initial field; 



3. solve the linearized equations to obtain an im- 
proved configuration; and 

4. iterate until the procedure converges. 

The third step above also involves a renormalization: the 
improved configuration after every Newton step is locally 
scaled to ensure that at every lattice point the constraint 
( ^.2[ ) is satisfied. This procedure speeds up the conver- 
gence to the final solution. The numerical solutions were 
found for m — 1 = g. 

The choice of initial configuration is guided by the 
known behavior of the periodic instanton solution near 
E ~ Q and E = Egpi^. Consideration of the zero mode 
in the vicinity of the sphaleron when /3 approaches /3_ 
indicates that we should look for non-trivial periodic so- 
lutions of the field equations in the neighborhood of the 
sphaleron by perturbing the known static solution along 
the zero mode direction, viz. 

ntspk)i^) ^ «(.p/.)W + S cos{\u;_\t)u'L{x) (5.11) 

where 5 is a small parameter which goes to zero at P — 
[3- . Since the classical action function on this solution is 
identical to Hamilton's principal function, the energy of 
the solution is given by 



dS{f3) 
dj3 



= E, 



(5.12) 



which must agree with the sphaleron energy Egph at 
f3 = (3-. For E just below Egp^ an initial configuration 
of the form (5.11) with (5 ~ 0.5 works well. The conver- 
gence of the Newton scheme is quadratic, and depending 
on the initial starting point, the desired accuracy (an er- 
ror tolerance of 10~^° or better) is typically reached in 
several iterations. Error tolerance limits were set to limit 
both the maximum violation of the equations of motion 
at a given lattice site as well as the global error obtained 
by summing the absolute value of, and then averaging, 
the errors on all lattice sites. Once the solution has been 
found at a given /3, we step down (or up) in [3 and use 
the previous solution as the initial trial configuration for 
Newton's method at the next value of (3. If the step size 
in /3 is not too large this procedure enables us to effi- 
ciently generate periodic instanton solutions over a finite 
range of (3. 

The solution of the linearized equations in step (3) 
of the algorithm can be performed by a direct inver- 
sion of the matrix of small fiuctuations about the trial 
configuration. (This was the method used in our earlier 
work on this problem and is described briefiy be- 
low.) Negative eigenvalues of this matrix are treated on 
the same footing as positive eigenvalues, and pose no spe- 
cial problem for Newton's method, which is important for 
the present application. Because of the boundary condi- 
tions which fix the translational and rotational symme- 
tries there are no zero eigenvalues of the matrix near the 
desired solution, which otherwise would be disastrous for 
this method. 
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For an I X J spacetime grid the matrix to be inverted 
has I^J^ elements. In general, it takes of order I'^J^ op- 
erations to invert such a matrix. However, the sparseness 
of the matrix makes it more efficient to follow a proce- 
dure of forward elimination and back substitution along 
either the a: or r directions instead. That is, starting with 
one edge of the region in Fig. 6, such as r = 0, we can 
solve for each r slice of the grid in terms of the successive 
two T slices until the edge r = /3/4 is reached, where the 
boundary condition determines the unknown quantities. 
Then we reverse direction and solve for the unknowns 
on the previous r slices successively. This allows for the 
matrix to be inverted in order IJ'^ operations (or I^J op- 
erations if the X direction is chosen) , and speeds up the 
algorithm considerably. In practice, a grid size of order 
100 X 100 can be handled on a typical workstation, which 
provides reasonable accuracy for configurations not too 
far from the sphaleron. However, this method rapidly 
becomes inefficient when large grids are required. 

The strategy adopted in this paper is to solve the lin- 
earized equations using direct matrix methods such as 
conjugate gradient. The present code allows for the use 
of several such solvers (all capable of dealing with non- 
symmetric, complex matrices). The results reported here 
were obtained using either the quasi-minimized version 
of Sonneveld's conjugate gradient squared (CGS) algo- 
rithm 1 22 1 or Van der Vorst's bi-conjugate gradient with 
stabilization (BICG_STAB) H]. The advantage of these 
methods over matrix inversion is a much less stringent 
memory requirement (down by a factor of / (or J) com- 
pared to matrix inversion) and while of the same order 
algorithmically, typically the prefactor is much smaller. 
Moreover, the efficient use of a parallel supercomputer 
becomes possible. 

Global grid refinement was used to speed up the 
method even further, as well as to effectively increase 
the order of the discretization error. The technique used 
was to begin with a certain grid size, find the periodic 
instanton, and then to double the grid size (keeping the 
physical volume constant) and use an interpolated ver- 
sion of the solution on the smaller grid as the initial guess 
solution for the bigger lattice. This procedure works as a 
sort of preconditioner for the matrix solver and since one 
has the results from the coarser lattice, Richardson ex- 
trapolation can be used to reduce the discretization 
error, effectively improving the second order finite differ- 
encing error to third order (Fig. 8). A good check of the 
error control achieved by the code is energy conservation. 
Energy error is maximum near r ~ /3/4 where the solu- 
tion is far from the vacuum and varying rapidly in both 
time and space. An example of the improvement in en- 
ergy error as a function of grid size (before Richardson 
extrapolation) is shown in Fig. 9. 

In practice, the number of (matrix solver) iterations 
hardly ever went beyond tens of thousands even for grids 
as big as of order 1000 x 1000 and acceptable accura- 
cies (error tolerances of order 10~^°) were reached in less 
than four doubling steps starting with configurations of 



order 100 x 100. When doing runs for different values of 
/?, the solution for the previous /S on the biggest grid was 
sampled to produce a guess solution for the new /? on 
the smallest grid. This allowed efficient scanning of the 
desired range of f3 values. As will be seen below the nu- 
merical results are extremely accurate all the way down 
to small values of /3 where the perturbative results of the 
previous section begin to apply. 




.2 2.3 2.4 2.5 2.5 2.7 2.B 2.9 3 



P 



FIG. 8. Action of the periodic instanton as a function of P 
illustrating the improvement due to extrapolation. The two 
crosses at each (3 are the numerical results from the largest 
and next-to-largest lattice at that /3 (top and bottom, respec- 
tively). The square is the extrapolated point. Accuracy of 
the extrapolation was checked by extrapolating up from a 
still coarser, i.e. next-to-next-to-largest lattice and compar- 
ing with the actual result at the largest lattice: there was very 
good agreement in all cases. 
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FIG. 9. Energy of the periodic instanton as a function of 
the time-slice at f3 — 2.3. A small range of r is shown to 
illustrate the improvement in energy conservation as the grid 
is refined (over four doublings of the grid in this case). Ex- 
trapolation improves the result one step further (not shown). 

Without any further optimizations, the new method 
is already efficient and accurate, yielding a real periodic 
instanton solution for energies in the range of .3 to 8 in 
units of m/g^. Smaller energies require finer lattices to 
obtain the same accuracy, since the instanton size be- 
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comes smaller rapidly with decreasing (3. Adaptive grid- 
ding using either nonuniform grids or selective refinement 
is possible within the method but was not used as a mean- 
ingful comparison with perturbation theory was already 
possible with the largest grid sizes that were used. (A 
selective refinement algorithm has been implemented for 
application to future, more demanding problems.) 

The results of the calculations are summarized below. 
They were obtained on grids as large as ^ 1200 x 2400 in 
spatial boxes with sizes varying from L = 8 to L = 16 in 
units of (for configurations nearer to the sphaleron 
smaller grids provided sufficient accuracy). For /?>/?_, 
as discussed in the previous section, a real solution does 
not exist and the solution becomes complex. This com- 
plex solution was also found without difficulty by our 
numerical method. Numerically, the crossover from the 
real to complex solution occurred between f3 = 3.62 (real) 
and P = 3.63 (complex) consistent with /3_ = 3.628 (Fig. 
4). The action of the complex solution rises steeply at 
/3 « 4 and it becomes difficult to track this numerical 
solution for larger /3 (possible in principle, but not worth 
the computer time). The numerical results for action and 
energy versus /3, and W as a function of E are given in 
Figs. 10. On the scale of the plots, the numerical error 
estimates are insignificant and are not shown. Compari- 
son with the perturbative results of the previous section 
shows very good agreement at low values of /3 as ex- 
pected. This is also iUustrated in Figs. 10 by plotting S, 
E, and W over a restricted range within which pertur- 
bation theory is supposed to become accurate. 

From Figs. 10 and 11, it is apparent that the peri- 
odic instanton solutions interpolate smoothly between 
the perturbative region and the sphaleron. In particu- 
lar, the energy monotonically increases with [3 for the 
entire range of /3, [0,/9_] in accordance with (1.1). 




FIG. 10. (a) The action of the periodic instanton as a 
function of [3. Agreement between the numerics and the per- 
turbation theory of the previous section is excellent in the 
expected range. 



Perturbation Theory 
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FIG. 10. (b) The same as Fig. 10(a) but with a more re- 
stricted range for /3 in order to better display the comparison 
with the low energy perturbative expansion (dashed curve) of 
Eqn. Oi. 




FIG. 11. (a) The energy of the periodic instanton as a 
funct ion o f (3. The dashed curve is the perturbative result, 
Eqn. ^^29. 



Perturbation Theory 
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FIG. 11. (b) The same as Fig. 11(a) but over a more 
restricted range of (3. 
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FIG. 12. (a) -W as a function of the energy E. The change 
of sign in W occurs at E = Esph = 8. 




1 1.5 

E 



FIG. 12. (b) The same as Fig. 12(a) but over a more 
restricted range of (3 in order to display the comparison with 
the perturbative result (dashed hne). 

The numerically evaluated periodic instanton configu- 
rations are displayed in Figs. 13 at = 2.8. These config- 
urations go over to the vacuum solution at the boundary 
X = — L/2 and are strongly non- vacuum near t = /3/4, 
X — 0. As /3 is lowered, the periodic instanton size 
goes rapidly to zero (exponentially in the perturbative 
regime). At lower values of /?, the configurations remain 
qualitatively similar, but become ever more localized as 
functions of time and space. 

The numerically computed periodic instantons can be 
compared to the pertur bativ e solution ( [3.21 ) with the in- 
stanton size p given by (4.27). Since direct comparison is 
unwieldy and difficult to visualize, we chose to compare 
the numerical and perturbative values for \w\ as a func- 
tion of r at a; = 0. This comparison at /? = 2.2 is shown 
in Fig. 14. As is to be expected, the overall agreement 
at this low value of /? is good. Evidently, perturbation 
theory underestimates the interactions in the dilute gas 
chain: the numerical solution has a longer tail and is less 
sharp than the perturbative configuration near r = /3/4. 





FIG. 13. A typical numerically obtained periodic instan- 
ton solution (here at /3 = 2.8, L = 12). Shown in the three 
figures are the Cartesian components of the unit vector , 
and n'^ as functions of r and a; on a 73 x 153 grid (r x x). As 
one can see, at the turning point r = the field configuration 
is non-vacuum, and the components of the vector n" cover 
the shaded patch of the sphere in Fig. 7. 
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FIG. 14. Comparison of the perturbative (solid line) and 
numerical results (dots) for \w\ as a function of r at a; = 0, 
/3 = 2.2. Only a small piece of the lattice near fH/A — .55 is 
shown for clarity. 




D^<f\-' + X[\<f\'-l-v 



(6.4) 

It is well-known that the instanton solution of this model 
is an Abrikosov-Nielsen-Olsen vortex of the form [ p5[ , 

1 

9 



$ = ve"*'H{r) . 



(6.5) 



In the special case that — 2Xv^ is taken equal to 



g V , the vortex may be found analytically, and 



has the action, Sq = irv'^ . For other values of the ratio 
Mh /Mw the vortex solution has to be found numerically 
p6| . The sphaleron and its negative frequency mode are 
both known analytically for arbitrary values of this ratio 



VI. SUMMARY AND DISCUSSION 

The main conclusion of our study of periodic instan- 
ton solutions in the quantum pendulum and the 0(3) 
sigma model is that two qualitatively and quantitatively 
different behaviors are possible for these solutions, with 
correspondingly different physical consequences for finite 
energy and temperature transitions. 

In the first case (I) of which the pendulum is the pro- 
totype, the period (3 decreases with increasing energy, 
Esph/T^ < 2So, the periodic instanton solutions con- 
tribute to the finite temperature winding number transi- 
tion rate, and the crossover from quantum tunneling to 
thermal activation is smooth: 



E 



and 



dE 



E. 



sph 



< 25-0 



case (I) 



.1) 



In addition to the periodic pendulum which we have dis- 
cussed here in some detail, the symmetric double well 
potential, 



(6.2) 



belongs in this class, as may be easily checked from 
Esph = V{0) = ^2/8: - |r'(0)|V27r - Lj/2V2n, 
and the kink action So = 2uj/3. Hence, 



Esph ttV2 4 
-jr- = ^-^^ < 2^0 = -uj , 



(6.3) 



and indeed the periodic instanton solutions with the ex- 
pected behavior are easily found analytically in this ex- 
ample as well. 

In two dimensions the most well-studied model pos- 
sessing both instanton and sphaleron solutions is the 
abelian Higgs' model, viz. 



sph 



Hence at least in the case = 2A we find 




Es 



87rV2 



3(VT7-f 1) 



< 2So = 2ttv^ 



(6.6) 



(6.7) 



and the two dimensional abelian Higgs' model also falls 
into case (I). The periodic instanton solutions have been 
found numericall y in this model for Mw = Mh , and 
they do satisfy (6T). In accordance with the intuition 
gained by our study of the sigma model, this behavior 
is the result of the absence of any conformal invariance 
or scale parameter akin to p in the instanton solutions 
of the abelian Higgs' model, and we would expect this 
model to behave in the same way for all finite values of 
the ratio Mh/Mw- The numerical results for the vortex 



action in the literature are consistent with (3/7) for all 
values of this ratio. 

The 0(3) sigma model is the only example studied so 
far in any detail which behaves differently, i.e. 



dE 
'did 



> 0; Esph/T_ > 2So . case (II) (6.8) 



In this second case (II), the periodic instanton solutions 
contribute to the finite energy (but not the finite temper- 
ature) winding number transition rate, and the crossover 
from quantum tunneling to thermal activation is sharp, 
taking place at Tcr > T- . As we have seen the key phys- 
ical difference between cases (I) and (II) appears to be 
conformal invariance and the scale parameter p of the 
instanton configurations. 

These same features are shared by the 0(3) sigma 
model and four dimensional gauge theories with sponta- 
neous (and conformal) symmetry breaking arising from 
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the Higgs' sector. Indeed, it is interesting to compare the 
numerical results for the sphaleron energy and negative 
frequency eigenvalue [ p9| , 



Mh 


Esph 




EsphOiW 
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± .O J-O 




0.1 


3.1384 
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16.18 


1.0 
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3.0 
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4.0 
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12.59 


5.0 


4.4375 


5.405 


11.99 


7.0 


4.6115 


7.176 


10.82 


10.0 


4.7805 


11.21 


8.97 



(6.9) 



to twice the zero energy instanton action 25*0 of the pure 
S'{/(2)-Higgs doublet model in four dimensions (which 
is the standard electroweak model with Oy/ — and 
no fermions). Since 2SQaw = 47r = 12.566... we 
observe from the last column of the table above that 
Esph/T- > 2So for all Mh/Mw < 4.0 but that the 
inequality is reversed for all Mw greater than a value 
slightly larger than AMw- Hence we conclude that the 
four dimensional S'J7(2)-Higgs' doublet model falls under 
case (II) for Mr < ^Mw- 

Based on our study of instanton solutions in this paper 
we are led to suspect that the four dimensional gauge the- 
ory is similar to the 0(3) model, for not too large Higgs' 
self-coupling. Indeed, since SU{2) pure gauge theory is 
conformally invariant, there is an instanton scale parame- 
ter p, which is driven to zero by the addition of the Higgs' 
sector with its conformal breaking expectation value v. 
Hence the instanton or anti-instanton solution no longer 
exists in isolation for non-zero Higgs' mass Mh = V2Xv. 
However, the interaction between instantons and anti- 
instantons is again attractive, so that at finite period- 
icity /3 there is an opposing interaction, and it is possi- 
ble to balance the two interactions at a definite value of 
p(/3), just as in the 0(3) model considered in this pa- 
per. As in the sigma model, in the spontaneously broken 
SU{2) gauge theory the period of the periodic instanton 
increases with energy at low energies, where perturbation 
theory is applicable 0. Thus we conclude that the four 
dimensional SU (2)-Higgs gauge theory falls into case (II) 
for not too strong Higgs' self-coupling, and has a sharp 
crossover at Tcr > T- from quantum tunneling to ther- 
mal activation in its winding number transition rate. 

Apparently the behavior changes at large Mh/Mw, as 
the model becomes more "Higgs-like" and the instanton 
solutions resemble the scale invariant SU (2) instantons 
less and less. At Higgs' self-coupling A > Snaw (approx- 
imately) the S'[/(2)-Higgs doublet model would appear 
to lie in case (I). At higher values of A additional static 
sphaleron solutions bifurcate from the simplest spheri- 
cally symmetric one, and the semiclassical approxima- 
tion which requires A <C 1 eventually ceases to be valid 
p9[ . These results have been obtained at zero U{1) mix- 
ing angle 9w- Although sphaleron solutions have been 



constructed at non-zero 9w |@ , the negative eigenmode 
and eigenvalue have not been given in the literature to 
our knowledge, and therefore we are not yet able to draw 
any firm quantitative conclusions in the case of non-zero 
6w- However, the behavior of the bosonic sector of the 
standard model with 9w near its physical value is not 
expected to be very different from the behavior at zero 
9w ■ In any case, we expect that the physical electroweak 
theory behaves like the 0(3) model for low to moder- 
ate Mh/Mw in case (II) but the situation with the ad- 
ditional Higgs coupling A in 3 + 1 dimensions is more 
complicated and deserving of detailed investigation. 

To conclude, the classical solution space of even rel- 
atively simple models such as the 0(3) model we have 
studied in this paper appears to be very rich, and the 
more intricate details of the periodic solutions to four di- 
mensional spontaneously broken Yang-Mills theories re- 
main to be investigated. Both the analytic and numerical 
approaches applied here to the 0(3) model will extend 
to the electroweak theory in a straightforward way, and 
we believe that this program is worth carrying through 
to completion. Besides the intrinsic interest of new solu- 
tions of the Euler-Lagrange equations of a physical gauge 
theory, this study, possibly extended into the complex 
domain appears to be the only viable method at our dis- 
posal to investigate anomalous baryon number violation 
at finite energy in the standard model. 
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